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Abstract — The problem of dirty paper coding (DPC) over 
the (multi-antenna) fading dirty paper channel (FDPC) 
Y = H(X + S) + Z is considered when there is imperfect 
knowledge of the channel state information H at the trans- 
mitter (CSIT). The case of FDPC with positive definite 
(p.d.) input covariance matrix was studied by the au- 
thors in a recent paper, and here the more general case 
of positive semi-definite (p.s.d.) input covariance is dealt 
with. Towards this end, the choice of auxiliary random 
variable is modified. The algorithms for determination of 
inflation factor proposed in the p.d. case are then general- 
ized to the case of p.s.d. input covariance. Subsequently, 
the largest DPC-achievable high-SNR (signal-to-noise ratio) 
scaling factor over the no-CSIT FDPC with p.s.d. input co- 
variance matrix is derived. This scaling factor is seen to be 
a non-trivial generalization of the one achieved for the p.d. 
case. Next, in the limit of low SNR, it is proved that the 
choice of all-zero inflation factor (thus treating interference 
as noise) is optimal in the 'ratio' sense, regardless of the co- 
variance matrix used. Further, in the p.d. covariance case, 
the inflation factor optimal at high SNR is obtained when 
the number of transmit antennas is greater than the num- 
ber of receive antennas, with the other case having been 
already considered in the earlier paper. Finally, the prob- 
lem of joint optimization of the input covariance matrix and 
the inflation factor is dealt with, and an iterative numerical 
algorithm is developed. 

Index Terms — Dirty Paper Coding, Inflation Factor. 



I. Introduction 

IN this paper, we continue with our study of the prob- 
lem of dirty paper coding (DPC) [T] over the fading dirty 
paper channel (FDPC) Y = H(X + S) + Z that was initi- 
ated recently in [2J. We consider the case in which there 
is imperfect and perfect knowledge of the fading process 
H at the transmitter (CSIT) and the receiver (CSIR), re- 
spectively. This problem is important from the point of 
view of application of DPC to the Gaussian multiple-input 
multiple-output (MIMO) broadcast channel (BC) with im- 
perfect CSIT. In [2J, we considered the case of input covari- 
ance matrix being positive definite (p.d.). Those results are 
restrictive in a sense that if a beamforming and DPC-based 
transmission scheme is to be used over the Gaussian MIMO 
BC (beamforming-based schemes are popularly being stud- 
ied for the BC) then the results of [2] alone will not suffice. 
With this motivation, this paper expands the scope to in- 
clude a more general case of positive semi-definite (p.s.d.) 
input covariance. 

The most crucial part of the problem of DPC over the 
imperfect-CSIT FDPC problem is to determine the op- 
timal inflation factor. Towards this end, for the case of 
p.d. input covariance, we proposed two iterative, numeri- 
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cal algorithms. These algorithms are seen to perform well. 
However, these algorithms can not be used directly when 
the input covariance matrix is p.s.d. This is because the 
achievable rate expression derived there is valid only when 
the covariance matrix is p.d. We thus slightly modify the 
choice of auxiliary random variable (RV) (see [3J, [I] for 
definition), derive the achievable-rate expression for the 
p.s.d. case, and then generalize the previously-developed 
algorithms to apply to the case of current interest. Next, in 
Section [PVl we derive the optimal high-SNR scaling factor 
achievable over the FDPC with p.s.d. input covariance. 
This scaling factor has several interesting characteristics 
that are not observed in the p.d. case. Further, an impor- 
tant result regarding the low-SNR behavior of the achiev- 
able rate is proved. Numerical results obtained using our 
algorithms are furnished in Section fVl 

Finally, we consider the problem of optimizing the 
achievable rate of the no-CSIT FDPC over the input co- 
variance matrix. Such a problem is well-studied for the 
MIMO channel [4], [5], [6]. The problem becomes more 
interesting in the case of FDPC because of the dependence 
of the inflation factor on the input covariance matrix. As 
a result, the covariance optimization must now be treated 
as a problem of joint optimization over the inflation factor 
and the input covariance. We develop an iterative algo- 
rithm for this joint optimization in Section [VII 

Notation: For a square matrix A, tr(vl), rank(A), 
and | .A | denote its trace, rank, and determinant, respec- 
tively. For any general matrix A, A*, A + , null (A) denote 
its complex-conjugate transpose, pseudo-inverse, and null 
space, respectively. I m is m x m identity matrix. Ex(-) 
denotes expectation over RV X. 

II. Channel Model 

The txr FDPC is defined via equation Y = H(X + 
S) + Z. Here, X ~ CA/*(0, Ex) is the transmitted signal 
with the power constraint of tr(Ex) < P; the interference 
S ~ CW(0, £s) is assumed to be known non-causally at the 
transmitter but not at the receiver; Z ~ CA/*(0,E.z) is the 
additive noise; and X, S, and Z are independent. The 
channel matrix H is assumed to be known perfectly at the 
receiver, whereas at the transmitter, only an estimate H 
is known. Assume \E Z \ > 0. Let tr(S s ) = Q, tr(£ z ) = N; 
and define SNR = 

Let the rank of Ex be m < t. Then we can write Ex = 
TT* for some ixm matrix T, and therefore, let X = TX' 
for some X' ~ CAf(0, I m ). We thus choose the auxiliary RV 
as U = X' + WS, as opposed to the choice U = X + WS 
made in the p.d. case. Here m x t matrix W is the inflation 
factor. Note that this choice of U is still a generalization 
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of Costa's (original) choice. 

Using the capacity formula of [7] , we obtain the achiev- 
able rate as given by 



(1) 



R = E H \og{\H(TT* + E s )iT + E z | }- 



muiE^log 



I m + WT, S W* (T* + WT, S )H* 
H(T + E S W*) h(tt* + y, s )h* + j:z 



Under this choice of U, the perfect-CSIT optimal infla- 
tion factor is given by W opt = t*H*(HTT*H* + Y, z )- l H. 

We define the no-interference upper-bound C as the rate 
achievable over the FDPC in absence of interference, i.e., 
when Q = 0. Thus if the covariance matrix of Ex is used 
over the FDPC, then C = E H log ^ z+ l "^ H ' 1 ■ 

III. Determination of Inflation Factor 

The minimization in equation ((T|) is the problem of de- 
termination of inflation factor. As noted in [2], it is a non- 
convex optimization problem and it seems intractable to 
find an analytical closed-form solution for the optimal infla- 
tion factor. We generalize below our previously-developed 
algorithms to apply to the problem at hand. 

Algorithm 1: In this algorithn0, we minimize the objec- 
tive function stepwise, i.e., at each step, we minimize over 
one row of W while treating all other rows as constants. 
Let us first consider the minimization over . Similar to 
the minimization problem considered in [2] , in this problem 
also, only the first row and the first column of the block- 
partitioned matrix in {2}, which we call M, depend on 

a B* 

W(i). Therefore, we repartition M as ' 



B D , where 

o = M (11)l B* = [W (1) E S (W)W (Tfo+W^Xsm, and 
D is (m + r — 1) x (m + r — 1) matrix remained after ex- 
cluding the first row and the first column from M. Now 
\M\ = \D\(a - B*D~ 1 B). Clearly, the problem reduces to 
the minimization of Eh log(a — B*D~ 1 B) over W(jy To- 
wards this end, observe that the expression (a — B* D~ 1 B) 
is a quadratic in Wny We now upper-bound the objective 
function through Jensen's Inequality by moving the expec- 
tation inside the logarithm, then the technique of 'complet- 
ing the square' can be used to minimize over Wm. 

If m = 1, the optimal W is given by equation ([2]) at 
the bottom of the next page, otherwise proceed as follows. 
" F G ~ 



Partition D as 



J K 



, where F is (m — 1) x (m — I) 



square matrix, K is r x r matrix, and G* — J. Then we 
can obtain W(i) as given by equation |pJ]Q The required 
expectations can be evaluated numerically. 

To minimize over the k th row of W, a similar idea can 
be used. Thus, an iterative algorithm can now be set up, 

1 Notations used in this algorithm: For any matrix A, A^ denotes 
the first row of it; A^ its first column while A^n^ is its (1, 1) element; 
Am and A^ denote entire matrix A except for its first row and 
except for its first column, respectively. 

2 Certain matrices inverted in equations J2]l and (O are singular 
if and only if Eg is p.s.d. This can be fixed by writing £5 = T2TJ , 
similar to what is done when Ex is p.s.d. 



as done in [2], in which one iteration consists of step- wise 
minimizations over all the rows of W and these iterations 
are repeated until a good choice is found. Importantly, this 
algorithm converges [2]. 

Algorithm 2: In this algorithm, we find a stationary 
point of the objective function, i.e., solve for Karush- 
Kuhn- Tucker conditions of the optimization. Thus we 
need to solve an equation -^^ H i H log \M\ = 0. It can 
be proved, as done in [2], that the above equation re- 
duces to Ejj.fiiAiW + A%H)E S = 0, where [A t A* 2 ]* = 
M~ x [I m 0]*. Without loss of generality (i.e., even if 
|Es| =0), we can consider a solution of the form W — 
-(E H ^ H A l y 1 E H ^(A^H) = g(W) for the above equa- 
tion^]. Using this solution, an iterative numerical algorithm 
can now be set up in which at the n th iteration, we set 
VF (n) = g(W (n ^ 1) ). 

IV. High and Low SNR Analysis 
A. High SNR Scaling Factor 

The scaling factor of min(t,r) is achievable over the no- 
CSIT FDPC, irrespective of Eg, when Ex is p.d. and it is 
equal to that of the no-interference upper-bound [2] . In the 
p.s.d. case, the high-SNR scaling factor may not necessar- 
ily be equal to that of the no-interference upper-bound and 
it also depends on Eg. Interestingly, such characteristics 
are never observed in the p.d. case. 

Theorem 1: Over the no (and hence, partial) CSIT 
FDPC with p.s.d. input covariance matrix, the 
largest DPC-achievabl^] scaling factor is given by 
min{r, rank(Ex + Eg)} — min{r, rank(Ex + Eg) — to}, pro- 
vided the ratio p is held constant as P — > 00, and the fad- 
ing process H is such that for any p.s.d. E, rank(i/Ei/*) = 
min(r,rank(E)) with probability 1. 

Proof: Let us reformulate the achievable-rate expres- 
sion as follows: define T via E x = PTT* so that X' ~ 
CN{Q,PI m ). Thus, we obtain R as given by equation 
([4} at the bottom of the next page where Eg = -PE' S . 
In order to maximize the achievable scaling factor, we 
need to choose W such that I) log \P(I m + WT,' S W*)\ 
does not achieve the high-SNR scaling factor larger than 
m, and 2) the term E H log |E Z + PHAH*\, where A = 
TT* + E^ - (T + E' s VF*)(/ m + WY,' s W*)- l {T* + WE' S ), 
scales minimally in the high-SNR regime. 

Let us analyze the second condition. It can be 
proved that A is p.s.d. for any choice of W . This 
and the assumption made about H imply that the 
term E H log |E Z + PH AH* \ would scale at high-SNR as 
min(r, rank(A)) logSNR. Therefore we need to minimize 
rank(j4) over W . To this end, we prove that rank(^4) > 
rank(TT* + Eg) - m. This lower-bound on rank(A) is 
achievable with the choice W = T + . Also, this choice of 
W satisfies the first condition. ■ 

3 Matrix A± is invertible with probability 1. 

4 The high-SNR scaling factor derived in Theorem 1 is optimal 
under the choice of DPC-based transmission strategy, but need not be 
the best achievable under all transmission strategies. This is because 
Costa's choice of auxiliary RV, which is considered here, is not known 
to be optimal under imperfect CSIT. 
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The assumption regarding the fading process H made 
above warrants a discussion. This assumption is fairly 
general and many fading distributions of practical inter- 
est would fall in this category. For example, if the proba- 
bility that null(iJ) contains certain fixed subspace (of the 
t-dimensional space) is zero then such a fading distribu- 
tion would fall in the category, and the commonly studied 
Rayleigh fading satisfies this condition. If this assump- 
tion is not made, the scaling factor would turn out to be 
dependent on the fading distribution. 

Let us now discuss the cases when the scaling factor of 
the no-CSIT FDPC and that of the no-interference upper- 
bound are equal. When t <r, two are equal (to m), irre- 
spective of the choices of Ex and Eg. However, unlike the 
case of p.d. covariance matrix [2], now, under p.s.d. in- 
put covariance, the difference between the no-interference 
upper-bound and the achievable rate may not necessarily 
go to zero in the limit of high SNR. Now consider t > r. 
If t > r > rank(Ex + Eg), two scaling factors are equal. 
On the other hand, when t > rank(Ex + E<?) > r, two are 
equal only if rank(Ex + Eg) = m. 

It is worth mentioning here that whenever r > m = 
rank(Ex + Eg), it can be shown that the difference be- 
tween the no-interference upper-bound and the achievable 
rate goes to zero in the limit as P — > 00. This follows from 
Theorem 2 of 0. 

Note that Proposition 3 in [2] about the D PC- achievable 
scaling factor of the no-CSIT Gaussian MIMO BC relies 
on the use of p.d. input covariance matrix for the last user. 
If a p.s.d. covariance matrix is used instead, we may get a 
scaling factor lesser than min(i,r). 

B. Low SNR Analysis 

Theorem 2: With the choice of W = 0, the ratio of 
the rate achievable over the no-CSIT FDPC to the no- 
interference upper-bound goes to 1 in the limit of low SNR 
if the ratio % is kept constant. 

Proof: With W = 0, we have 
R = Eh log \T, z + H(^ x + ^s)H*\-E H log |E Z + H^ S H* | ; 
and the no-interference upper-bound as 



would 

be of the form (R,C) = a(R,c)P + /3 (fliC) P 2 + o(P 2 ), where 

Hence, if 



a 



(R,C) 



d(R,C) 
dP 



:,c) 

d 2 (R,C) 
dP 2 



P=0 



C = E H log |E Z + H^ X H* I - log |£ z 

Taylor-series expansions of 7? and C about P 

(R,C) 

, and /3(k,c) 

= otc, we must have limp^o % = 1- 
Let us first derive olr. Let Ex = PE^-, Eg = 
[Ez]y = Zij, [HT,' x H*]ij — aij, and [HT,' s H*]ij 
Then for M = E z + PH(E X + E' S )H*, we have 

— log \M\ = - — - V I M W I where 



PE'c 



5- 



[M^] jk =z jk +P(a jk + b jk ) 
= {ajk+bjk) 



if j 7^ i, and 
if j = i. 



i th row re- 



Here, is nothing but the matrix M with 

placed by its derivative with respect to P. This step can be 
verified from the so-called 'ugly' formula of determinants. 
Letting N = ~E Z + PPE S P*, we get a R = 



where [Q {l) )jk = z jk 
= a.jk 



if i ^ j, and 
if i = j. 



This follows because at P — 0, all rows of AfW and those of 
are equal except the i row. The result then follows 
by noting that a,j k = [HT,' x H*]j k . I 
Thus, in the limit of low SNR, it is optimal, in a 'ratio' 
sense, to treat the interference as noise. Note that this re- 
sult holds irrespective of whether p.d. or p.s.d. covariance 
matrix is used. 

C. High-SNR Analysis: p.d. covariance matrix 

We include in this subsection one important result re- 
garding the p.d. case. It is proved in [5] that the choice 
W — It is optimal in the high-SNR regime if t < r. We 
now prove a similar result for the remaining case of t > r. 



II 7 ' 11:1 I ' + E S )P* + E z ]" 1 7J}E s (Es - E S E{PT*(P(TT* + Es)P* + 'E Z )~ 1 H}E S )~ 1 , 



(2) 



W (1) = ({I m ) {1) T*E(H*J)W (1) -£s + (I m )(i)T*E(H*KH)i:s) x (E s - E S (W (I) )*E(P)W (I) E S 
-E S E(P* J)W a) E s ~ E s (PF ( i))*E(GP)E s - E S E(P* A-P)E S ) _1 , • • • where E(-) = E ff|A (-). (3) 



R = E H log 



\PI m \ |E Z + PH(TT* + E' 5 )P* 



\P(I m + VFE' S M/*)| |E Z + PH{TT* + E' s - (T + H' s W*)(I m + W'E' S W*)~ 1 (T* + VFE^)}P* 



(4) 



R = E H log- 



■H(H 



x 



E«)P* 



E x + WT, S W*\\T, Z + PH{Z' X + E' s - (E x + Y>' S W*){Y>' X + WE'sW*)- 1 (£' x + WY,' S )}H* 



(5) 



-1 



Theorem 3: For the no-CSIT FDPC with t > r and p.d. 
input covariance matrix, the inflation factor optima^ in 
the high-SNR regime is given by W — I t + AW, where 
AW is any matrix such that AWT,s — 0, provided the 
ratio f is kept constant, and H is such that for any p.s.d. 
S, rank(iJ£iJ*) = min(r,rank(£)) with probability > 0. 

Proof: The achievable-rate expression from [2] can be 
rewritten to get equation ([5]) at the bottom of the previ- 
ous page. Let A(W) = T,' x + Z' S - (E^ + Z' S W*)(E' X + 
WH'sW*)- 1 ^ +WZ' S ). The idea is to characterize the 
set of all choices of W that achieve the optimal scaling. 
It turns out that the choices of W described in the state- 
ment of the theorem are the only choices that achieve the 
optimal scaling, and hence follows the theorem. 

5 Since we considering the optimality of W, let Wp be describe 
the inflation factor used at power P. We assume in this theorem that 
limp^^ Wp exists. 



Now, R scales optimally if and only if (iff) W is such 
that 1) the factor log |£x + WT,$W* \ does not achieve the 
high-SNR scaling factor larger than t, and 2) the term 
E ff log|E z + PHA(W)H*\ does not scale in the high-SNR 
regime. To show the dependence of A on the transmitted 
power P, let us write A(W) as A(W P ). Note that A(W P ) 
is p.s.d. for any choice of Wp. 

Let us analyze the second condition. The assumption 
made regarding H implies that the second condition holds 
iff limp_ >00 tr{A(Wp)} = 0. Since the matrix inverse is a 
continuous function (note (T,' x + WpT,' s Wp) is always p.d.) 
and we have assumed that the limit of Wp exists, the above 
implies the existence of W = W^ — limp^oo Wp such that 
tr{A(W)} = tr{A(W r 00 )} = 0. Since A(W) is p.s.d. for any 
W, this implies ^4(M^) = 0. Therefore, let us characterize 
the set of all W's for which y^VF) can be zero. 





Fig. 2 

Achievable Rate vs. SNR for 3 x 2 System, (a) rank(Sx + £s) > rank(Ex)- (b) rank(Sx + £s) = rank(Sx)- 
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At optimal W, for which A(W) = 0, the matrix (E' x + 
WE' S ) must be invertible since matrices (E' x + T,' s ) and 
(Y,' x +WY,' S W*) are p.d. Therefore at optimal W, E x + 



WY,' S W* - (E' x + WZ' S )(Z> X + V S )-\X' X 



E' S VF*) = 0. 



The left hand side of this equation is a quadratic in W . 
Hence, it can be written as (Wa — b)(Wa — b)* + c for 
some appropriately chosen a, 6, and c. Equating terms, 
we get aa* = Eg — Eg(Ex + an d ab* — Eg(Ex + 

Es) _1 Ex- Note that aa* = ab* and hence without loss of 
generality, we can set a = b. This implies c = 0. 
Case I: Eg is p.d. aa* is p.d. Then there exists a unique 
solution W = It- 

Case II: Eg is p.s.d. aa* is p.s.d. In this case, 

W — h is one of the possibly many solutions. Let the 
optimal solution be W = I t + AW. Then, we must have 
(AW)aa* (AW)* = 0. But rank(aa*) = rank(Eg), and 
hence null(aa*) = null(Eg). This implies AWE S = 0. ■ 
Since W always gets multiplied by Eg in the rate ex- 
pression, we can set AW = 0, without loss of generality. 

V. Numerical Results 

In Figs. [H to [31 we consider real-valued FDPCs. Here, 
the elements of H are i.i.d. ^A/"(0,1). We quantize each el- 
ement of H separately using a simple 'equally spaced level' 
quantizer (quantization bins are of equal length except for 
the first and the last bins which extend to infinity) and the 
spacing is determined using data from [5] . B denotes num- 
ber of feedback bits per element of H . 'no-S-ub' denotes 
the no-interference upper-bound. Also in these results, we 
use Algorithm 1. Figs. [1] and [2] have two parts. The only 
difference between the FDPCs of part (a) and (b) is that 
the choice of Eg is varied 0. 

In Figs. Q]and[2l we consider 2x2 and 3x2 FDPCs, re- 
spectively. In both figures, we see a significant increase in 
the achievable rate with increasing feedback. Observe that 
for the 2 x 2 FDPC, the high-SNR scaling factor of R and 
that of the 'no-S-ub' are equal whereas for the 3x2 FDPC 
of Fig. [21 the scaling factor of R is dependent on the choice 
of Eg. Interestingly, in Figs, [lib) and [2(h), the difference 
between the no-interference upper-bound and the achiev- 
able rate can be seen to tend to zero with increasing SNR. 
The high-SNR behavior of the achievable rate observed in 
these figures is in accordance with the analytical results 
proved earlier. 

In Fig. [31 we consider the low-SNR behavior. At suf- 
ficiently low SNR, the choice W = can be said to be 
near optimal. This fact is formulated mathematically and 
proved in Theorem 2. 

In Fig. [4l we present comparison of two algorithms. We 
take complex-valued systems with H Rayleigh faded and 
also consider the case of no CSIT. It is observed rather 
surprisingly that these algorithms typically yield W's that 
achieve almost equal rate and therefore, we present here 
some cases of difference. For FDPC 1, m = 1 and Algo- 
rithm 1 uses the closed-form solution obtained in equation 
©. Curves for FDPC 1 suggest that Algorithm 2 yields a 

6 All parameters necessary to generate the plots enclosed in this 
paper can be obtained by e-mailing the first author. 




Fig. 3 

Low-SNR Analysis: Achievable Rates vs. SNR. 



solution that is slightly better than this closed-form solu- 
tion of @. The main advantage of Algorithm 1 over the 
second is that it is faster in terms of number of iterations 
required. 

VI. Joint Optimization over the Input 
Covariance Matrix and the Inflation Factor 

We solve this problem of joint optimization under an 
additional constraint on rank(Ex). 

Since using Algorithms 1 and 2, we can determine W for 
a given choice of Ex, let us first consider optimizing over 
Ex for a given W. Let the rank of Ex be constrained by m. 
Under this condition, we write Ex = TT* for some t x m 
matrix T, and therefore, rank(Ex) < m. This also ensures 
that Ex is p.s.d. Thus we have the following optimiza- 
tion problem: maxp i?, subject to the power constraint 
tr(TT*) < P. 

To derive the necessary conditions for the optimization, 
we form the lagrangian J and set ^ = 0. We have 



J = Eh log 



\Z z + H{TT* + Y<s)H* 



I m + WT, S W* (T* + WT, S )H* 
H(T + E S W*) Ez + H(TT* + E5)iJ" 

- Atr(TT*), 



where A is the lagrange multiplier. 

As mentioned before, the problem of covariance opti- 
mization for the MIMO channel has been studied in 0] , [S] , 
[B]. In dj, a suboptimal solution is provided; whereas in 
[S] and [BJ, the authors solve for the Karush-Kuhn- Tucker 
conditions, and to this end, differentiate the lagrangian 
with respect to each element of T separately. This tech- 
nique would be very cumbersome for the problem at hand. 
The key here is to differentiate the lagrangian with respect 
to matrix T directly using the techniques of matrix dif- 
ferentiation. Due to lack of space, we omit the details of 
finding derivatives and directly state below the necessary 
conditions for the above optimization: 



XT = E H \ FN^HT + [0 H*]D; 



■* m 

HT 



g{T,w). 
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Fig. 4 

Comparison of Algorithms: Achievable Rates vs. SNR. 



This fixed-point equation allows us to set up an iterative 
algorithm for the joint optimization. 

1. Start with some initial choice for T, say T(°) . 

2. At the n th iteration: 

. For the choice T^ -1 ) of T, find W. 

• For VF" obtained above, find optimal T. Towards this 
end, we use the above equation, and set T^ n > — 
i(7(T(™ -1 ), W), where the required expectations are 
computed numerically and the lagrange multiplier 
A is evaluated so as to meet the power constraint. 

3. Repeat the above step until the increase in the achiev- 
able rate is negligible. 

Numerical Results: In Fig. El FDPC is the rate achiev- 
able over the FDPC with the jointly-optimized input co- 
variance and the inflation factor and similarly for curves 



labeled MIMO. FDPCib is the rate achievable using the 
scaled-version of identity as the input covariance matrix 
and the optimal inflation factor (determined using Algo- 
rithm 1). If £ is the ratio of the largest to the smallest 
eigen- value of the input covariance matrix (obtained using 
the algorithm of joint optimization), then in Fig. 0(b) , 
'Eig. Ratio' denotes ^log 10 (x). 

In Fig. EJa), we consider correlated Rayleigh fading with 
separable correlations, whereas in Fig. E[b), we take ele- 
ments of H to be i.i.d. ~ Unif [0, 1] + jUnif [0, 1] where 
Unif[0,l] denotes uniform distribution over [0,1]. In these 
figures, the improvement due to spatial water-filling is ev- 
ident. In Fig. [5](c), we consider a rank-constrained opti- 
mization and H is Rayleigh faded. One can see the suf- 
ficiency of reduced-rank signalling at low SNR (rank 1 at 
— 10 dB, 2 at —5 dB) while the full-rank signalling performs 
significantly better (than the reduced-rank signalling) at 
high SNR of 30 dB. 



VII. Conclusion 

The determination of the inflation factor under partial 
CSIT was the major obstacle in the problem of apply- 
ing DPC to the partial-CSIT Gaussian MIMO BC. The 
results of this and our earlier paper [2] make significant 
progress to address this problem. At this time, the prob- 
lem of jointly optimizing the sum-rate of the partial-CSIT 
Gaussian MIMO BC over the covariance matrices of all 
the users (and the inflation factors) seems to be the next 
logical step in developing an algorithmic solution for com- 
puting good lower bounds on the achievable sum-rate of 
the Gaussian MIMO BC. 




Fig. 5 

Covariance Optimization, (a) Ratio of Various Rate Pairs for 3x3 System, (b) Ratios for 2x2 System, (c) Rank 

Constrained Optimization for 3x2 System. 
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